Metal-Insulator-Like Behavior in Semimetallic Bismuth and Graphite 
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When high quality bismuth or graphite crystals are placed in a magnetic field directed along the 
c-axis (trigonal axis for bismuth) and the temperature is lowered, the resistance increases as it does 
in an insulator but then saturates. We show that the combination of unusual features specific to 
semimetals, i.e., low carrier density, small effective mass, high purity, and an equal number of elec- 
trons and holes (compensation), gives rise to a unique ordering and spacing of three characteristic 
energy scales, which not only is specific to semimetals but which concomitantly provides a wide win- 
dow for the observation of apparent field induced metal-insulator behavior. Using magnetotransport 
and Hall measurements, the details of this unusual behavior are captured with a conventional multi- 
band model, thus confirming the occupation by semimetals of a unique niche between conventional 
metals and semiconductors. 

PACS numbers: 71.30.+h, 72.20.My, 71.27. +a, 73.43.Qt 



Elemental semimetals, such as bismuth and graphite, 
are intriguing materials to study because of their high 
magnetoresistance, low carrier density n, and high pu- 
rity. Due to small values of n, magnetic fields on the 
order of 10 T are sufficient to drive these semimetals 
into the ultraquantum regime, where only the lowest Lan- 
dau level remains occupied. In addition, light cyclotron 
masses, m* , for any field orientation in Bi and along 
the c-axis in graphite, result in higher cyclotron frequen- 
cies, lo c = eB/m*, which ensure that quantum magneto- 
oscillations can be observed at moderate temperatures. 
High purity facilitates the observation of well-resolved 
oscillation patterns. These features have made bismuth 
and graphite perhaps the two most popular materials for 
studies of quantum magnetic-field effects Q, Q . 

Recently, interest in magnetotransport in graphite has 
been renewed due to observations (similar to those shown 
in Fig. 1) of an effect that looks like a magnetic-field- 
induced metal- insulator transition 0, The metal- 
lic T-dependence of the in-plane resistivity in zero field 
turns into an insulating- like one when a magnetic field 
on the order of 10 mT is applied normal to the basal 
(ab) plane. Increasing the field to about 1 T produces a 
re-entrance of the metallic behavior [j| . It has been pro- 
posed that the low-field effect is due to a magnetic-field- 
induced excitonic insulator transition of Dirac fermions 
U |(J, whereas the high-field one is a manifestation of 
field-induced superconductivity . It has been also sug- 
gested that the apparent metal-insulator transition 
in graphite is similar to that in 2D heterostructures 
(although the latter is driven by a field parallel to the 
conducting plane). 

Similar metal-insulating like behavior is also observed 
in 99.9995% pure bulk bismuth crystals as shown in 
Fig. 2, where the resistivity is plotted as a function of 
temperature in successively higher magnetic fields. The 
crossover from "metallic" to "insulating" behavior has 
the same qualitative behavior in both semimetals. On 
lowering the temperature the resistance increases, as it 



does in a insulator, but then saturates towards field- 
dependent constant values at the lowest temperatures. 
These similarities invite an interpretation that ascribes 
this interesting behavior to properties shared by both 
graphite and bismuth, namely low carrier density, high 
purity, and an equal number of electrons and holes (com- 
pensation), rather than to specific properties of graphite: 
almost 2D nature of transport and a Dirac-like spectrum, 
as suggested in Refs. 0, 0> H, E| • 

In this letter we demonstrate this connection by ex- 
amining the magnitude and ordering of three character- 
istic energy scales: namely, the width of the energy lev- 
els Ti/t where r is the electron-phonon scattering time, 
the cyclotron energy Juu c , and the thermal energy fc^T. 
We provide theoretical justification for, and experimental 
confirmation of, the existence of a wide interval of tem- 
peratures and magnetic fields, defined by the condition, 



(1) 



In this interval, (a) the magnetoresistance is large, (b) the 
scattering rate is linear in T, and (c) Shubnikov de Haas 
(SdH) oscillations are not resolved due to the thermal 
smearing of Landau levels. We argue that the unusual be- 
havior of bismuth and graphite is due to the existence of 
a region, specified by the inequalities of Eq. Q), and also 
due to compensation between electron and hole charge 
carriers. Our experimental confirmation is centered on 
a detailed study of graphite in which we use the con- 
ventional theory of multi-band magnetotransport to 
extract the field- independent carrier density, n(T), and 
scattering time, t(T), from simultaneous fitting of the 
temperature and field-dependent longitudinal resistivity 
Pxx{T,B) (magnetoresistance) and transverse resistivity 
p X y(T, B) (Hall effect). We then show from this analysis 
that the inequality QJ, which is unique to semimetals, is 
satisfied over a broad temperature and field range. 

To illustrate the uniqueness of low-carrier-density 
semimetals, we compare them with conventional, high- 
density, uncompensated metals. To begin with, if the 
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Fermi surface is isotropic, a metal exhibits no magne- 
toresistance because the Lorentz force does not have a 
component along the electric current In anisotropic 
metals the magnetoresistance is finite and proportional 
to (lo c t) 2 in weak magnetic fields, i.e., for u> c t < 1. In 
stronger fields (uj c t 3> 1), classical magnetoresistance 
saturates, if the Fermi surface is closed |9]. In contrast, 
transverse magnetoresistance of a semimetalgrows as B 2 
both in the weak- and strong- field regimes j2( • 

The magnetoresistance [p(B) — p(0)]/p(0) is much 
larger in semimetals than in conventional metals. In 
addition to the saturation effect, described above, an- 
other important factor that limits the magnetoresistance 
in conventional metals is the higher scattering rates and 
thus smaller values of the u> c t product. The impurity 
scattering rate in semimetals is smaller than in con- 
ventional metals simply because semimetals are typi- 
cally much cleaner materials. The lower carrier den- 
sity of semimetals also reduces the rates of electron- 
phonon scattering in semimetals compared to that of 
conventional metals. For temperatures above the trans- 
port Debye temperature, which separates the regions of 
the T- and T 5 -laws in the resistivity, Q P D — 2hkps/kB, 
where kp is the Fermi wavevector and s is the speed 
of sound (both properly averaged over the Fermi sur- 
face), one can estimate the electron-phonon scattering 
rate as r _1 ~ (fcp-ao) (m* /mo) ksT /h, where ag is the 
atomic lattice constant, and m* and too are respectively 
the effective and bare electron masses jjjj- I n a con ~ 
ventional metal, kpa-o — 1 and m* ~ too. In this case, 
Q P D is of order of the thermodynamic Debye tempera- 
ture hs/ksao — few 100 K and t" 1 ~ kBT/Ti. In 
low-carricr-dcnsity materials (kpao <C 1), which typi- 
cally also have light carriers (m* -C too), Q p d is much 
smaller and thus r _1 <C ksT/Ti. Therefore, in a low- 
carrier-density semimetal there exists a wide interval of 
temperatures and magnetic fields, defined by the inequal- 
ities In contrast, there is no wide interval between 
Ti/t and ksT in a conventional metal [llj . 

An additional feature, crucial for interpreting the ex- 
perimental data, is that the Fermi energies of graphite 
(E F = 22 meV) 12| and bismuth (E F = 30 meV 
[holes] )0 are relatively low and the temperature de- 
pendence of the resistivity comes from two temperature- 
dependent quantities: n(T) and t(T). Purity of materi- 
als ensures that electron-phonon scattering is a dominant 
mechanism for resistance over a wide temperature range. 

Standard 4-probe measurements were carried out on 
single-crystal highly oriented pyrolytic graphite (HOPG) 
sample with a 2° mosaic spread, as determined by X- 
rays. The resistivity was measured using an ac (17 Hz) 
resistance bridge over the temperature range 5K-350K. 
In all the measurements, the magnetic fields were ap- 
plied perpendicular to the sample basal planes. Both 
p xx and p xy (Fig. 3) were measured in magnetic fields 
up to 1 T, although the analysis (solid lines) was limited 
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FIG. 1: Temperature dependence of the ab-plane resistivity 
pxx for a graphite crystal at the c-axis magnetic fields indi- 
cated in the legend. The solid lines are the fits to the data 
using the six parameters derived from the three-band analysis 
described in the text. The shadowed region in the inset and 
its mapping onto the data in the main panel indicates the 
interval defined by Eq.Q. 



to B < 200 mT. A small field-symmetric component due 
to misaligned electrodes was subtracted from the p xy (B) 
data. 

We used a standard multi-band model Q to fit the 
data. Each band has two parameters: resistivity pi and 
Hall coefficient Ri = l/(fonj, where qi — ±e is the charge 
of the carrier. In agreement with earlier studies, we fix 
the number of bands to three . Two of these are the 
majority electron and hole bands, and the third is the 
minority hole band [2j. Although the third band is not 
essential for a qualitative understanding of the data, it is 
necessary for explaining the low-field fine features in p xy 
shown in the inset of Fig 3. The minority band makes 
a negligible contribution at higher fields due to its low 
carrier concentration. 

We fit p xx and p xy simultaneously by adjusting the 
six parameters independently, until differences between 
the fitting curves and the experimental data are min- 
imized. Because the majority carriers in graphite 
derive from Fermi surfaces that have six-fold rota- 
tional symmetry about the c-axis, we only need to use 
the 2x2 magneto-conductivity tensor a % with elements 
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to* /riie 2 Ti. The total 



, where pi 

conductivity, a — J2t=i ^ s s i m ply a sum of the con- 
tributions from all the bands and the total resistance is 

Qualitatively, the unusual temperature dependence 



3 




50 100 150 200 250 300 
T(K) 



FIG. 2: Temperature dependence of the longitudinal resis- 
tivity for a bismuth crystal at magnetic fields ranging from 
to f T in 0.1 T steps from top to bottom. The magnetic 
field and the current are along the trigonal and binary axes, 
correspondingly. In zero field, the resistance is approximately 
linear in temperature. 



displayed in Fig. 1 can be understood for a simple two- 
band case where p xx reduces to 



Pip2 (pi + Pi) + (p\R\ + p 2 Rj) B 2 
(pi + Pif + (Ri + R2) 2 B 2 



(2) 



Assuming that p\^ oc T a with a > 0, we find that for 
perfect compensation, i?i = — i?2 = \R\, Eq. © can be 
decomposed into two contributions: a field-independent 
term oc T a and a field-dependent term oc R 2 {T)B 2 /T a . 
At high T, the first term dominates and metallic behavior 
ensues. At low T, R(T) oc l/n(T) saturates and the 
second term dominates, giving "insulating" behavior. 

The actual situation is somewhat more complicated 
due to the T-dependence of the carrier concentration, 
the presence of the third band, and imperfect compen- 
sation between the majority bands. Results for the 
temperature-dependent fitting parameters are shown in 
Fig. 4 where band 1 corresponds to majority holes, 
band 2 to majority electrons and band 3 to minority 
holes. The insulating-like behavior of the carrier den- 
sity with a tendency towards saturation at low temper- 
atures is well reproduced. For the majority bands, 1 
and 2, the carrier concentrations are approximately equal 
and similar in magnitude to literature values |2j. The 
slope of the linear-in-T part of T^ 1 =a cxv ksT '/h with 
ctoxp = 0.065(3) (dashed line in Fig. 4, top panel) is con- 
sistent with the electron-phonon mechanism of scatter- 
ing. To see this, we adopt a simple model in which carri- 
ers occupying the ellipsoidal Fermi surface with param- 
eters m a b (equal to 0.055rno and 0.040toq for electrons 
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FIG. 3: Longitudinal resistivity p xx (top panel) and trans- 
verse resistivity p xy (bottom panel) versus applied field B for 
graphite at the temperatures indicated in the legend. The 
solid lines are determined by a fitting procedure that simul- 
taneously includes both sets of data into a three-band model 
described in the text. The inset in the bottom panel, with 
the same units on each axis, magnifies the low-field region, 
where the contribution from the third band with a lower car- 
rier density makes a major contribution that cannot be fit by 
solely using the two majority bands. 



and majority holes, correspondingly) and m c (equal to 
3mo and 6mo, correspondingly) interact with longitudi- 
nal phonons via a deformation potential, characterized 
by the coupling constant D (equal to 27.9 eV) . In this 
model, the slope in the linear-in-T dependence of r _1 is 
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given 
where m 

and holes, po=2.27 g/cm 3 is the mass density of graphite 
and s a b = 2 x 10 6 cm/s is the speed of sound in the ab- 



{m*fE F D 2 /p s 2 ab h\ 
0.21too both for electrons 
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FIG. 4: Temperature dependence of the fitting parameters 
(graphite) for the bands indicated in the legends of each panel. 
The near equality of the carrier densities of the two majority 
bands (lower panel) indicates good compensation at low fields 
(B < 200 mT) and the linear dependence of the scattering 
rate on T (upper panel, dashed line) with a slope of 0.065(3) 
in units of fcs/Ti is consistent (see text) with electron-phonon 
scattering. 



plane. (The numerical values of all parameters are taken 
from the standard reference on graphite @.) With the 
above choice of parameters, atheor = 0.052 for both types 
of carriers. This value is within 20% of the value found 
experimentally. Given the simplicity of the model and 
the uncertainty in many material parameters, especially 
in the value of D, such an agreement between theory and 
experiment is quite satisfactory. 

The solid lines through the data points in Fig. 1 are cal- 
culated from the temperature-dependent fitting parame- 
ters derived from our three-band analysis and plotted in 
Fig. 4. The shaded region (II) depicted in the inset of 
Fig. 1 represents those temperatures and fields that sat- 
isfy the inequalities in (1). In region (I) SdH oscillations 
can be seen at sufficiently low T. Our sample has a Dingle 
temperature of 3.0 K. The boundary between (I) and (II) 
reflects the rightmost inequality of (1) and is determined 
by the relation T > TieB/m*. The boundary between 
(II) and (III) reflects the leftmost inequality of (1) and is 
determined by the relation B > m* /er(T) where l/r(T) 
is obtained from experimental fitting parameters (Fig. 4). 
In the main panel of Fig. 1 we superimpose region (II) , 
again as a shaded area, on the p xx (T,B) data. Below 
the lower boundary oj c t < 1, and the magnetoresistance 
is relatively small. The upper boundary is determined 
by the locus of (B, T) points satisfying the rightmost in- 



equality of (1). Clearly region (II), constrained by the in- 
equalities of (1), overlaps well with the metal-insulating 
like behavior of graphite. Since the majority bands of bis- 
muth comprise three non-coplanar electron ellipsoids and 
one hole ellipsoid, a similar analysis for bismuth is more 
complicated and would require more space than available 
here. 

We thus conclude that the semimetals graphite and, by 
implication, bismuth share the common features of high 
purity, low carrier density, small effective mass and near 
perfect compensation and accordingly obey the unique 
energy scale constraints that allow pronounced metal 
insulating behavior accompanied by anomalously high 
magnetoresistance. At magnetic fields higher than dis- 
cussed in this paper (B > 1 T) we believe that the multi- 
band model is still appropriate and may provide an al- 
ternative explanation for the reentrant behavior observed 
by us and others |5j. 

Subsequent to the completion of this study, we were 
informed of recent work [T. Tokumoto et al, Solid State 
Commun. 129, 599 (2004)] which used a two-band model 
to explain the unusual behavior of p xx (T, B) in graphite. 
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